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Abstract 

Photoassociation of a pair of cooled atoms by excitation with a short chirped laser pulse creates a 
dynamical hole in the initial continuum wavefunction. This hole is manifested by a void in the pair 
wavefunction and a momentum kick. Photoassociation into loosely bound levels of the external well 
in Cs2 0^(6S+6P3/2) is considered as a case study. After the pulse, the free evolution of the ground 
triplet state a'^S^ wavepacket is analyzed. Due to a negative momentum kick, motion to small 
distances is manifested and a compression effect is pointed out, markedly increasing the density 
of atom pairs at short distance. A consequence of the hole is the redistribution of the vibrational 
population in the a^S^ state, with population of the last bound level and creation of pairs of hot 
atoms. The physical interpretation makes use of the time dependence of the probability current and 
population on each channel to understand the role of the parameters of the photoassociation pulse. 
By varying such parameters, optimization of the compression effect in the ground state wavepacket 
is demonstrated. Due to an increase of the short range density probability by more than two orders 
of magnitude, we predict important photoassociation rates into deeply bound levels of the excited 
state by a second pulse, red-detuned relative to the first one and conveniently delayed. 
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I. INTRODUCTION 



One of the main outcomes of excitation of a molecular system with a short laser pulse is 
a "hole" in the wavefunction of the initial ground state [2|- The impulsivee limit has been 
extended to the problem of photoassociation of ultracold atoms Ref. 3j. Due to the ultralow 
collision energy this limit can be applied to a well-defined initial stationary collision state. 
The emphasis of that study was on optimizing the photoassociation process and the creation 
of stable ultracold molecules. We should complete the story and follow the dynamics of the 
"hole" carved out of the initial continuum wavefunction. We have already identified the 
signatures of this " hole" : 

• 1) a void in the pair wave function of the initial state 

• 2) a momentum kick, initiating a motion of the hole 

• 3) creation of stable molecules , via population of the last bound levels of the ground 
electronic state. 

• 4) creation of pairs of hot atoms, via population of continuum levels with higher energy 

The aim of the present paper is to analyze in detail the hole created in the ground state 
continuum wave function, during the photoassociation pulse and after it. In Section [Tll we 
display results of numerical calculations on photoassociation of cesium with chirped laser 
pulses, showing the time-evolution of the wave packet in the initial state. We identify the 
dynamical hole, its motion, and the compression of the initial wavefunction, increasing the 
probability density at short internuclear distances. Possible applications of this compression 
effect to photoassociation with a second pulse are suggested. In Section [Till we present 
theoretical tools for the analysis of a dynamical hole through the expectation value of the 
current density vector. In Section [IV] we make use of such tools to analyze the numerical 
results of Section [ITl and their sensitivity to the pulse parameters (energy, duration, central 
frequency, sign a nd value of the chirp parameter). In Section |V] we discuss how to control the 
formation of a dynamical hole and the compression effect, and to optimize photoassociation 
into vibrational levels with low v values. 
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II. PHENOMENOLOGICAL DESCRIPTION OF THE DYNAMICAL HOLE AND 
OF THE COMPRESSION EFFECT: NUMERICAL CALCULATIONS IN CASE OF 
CESIUM PHOTOASSOCIATION. 



A. The physical system 

The photoassociation reaction 

2Cs(65, F = 4) + h{u;{t)) ^ Cs2(0^ (65 + 6P3/2); t^, J) (1) 

is employed as a primary example. In this reaction two ultracold cesium atoms absorb 
a photon red detuned from the resonance line to form a bound level {v, J) in the outer 
well of the excited potential curve 0^(65' + 6P3/2). The potential curves have been described 
elsewhere [s] and are displayed in Fig. [TJ We shall use vibrational numbering for the double- 
well potential : the external well supports at least 227 levels, from f = 25 up to = 256, 
the levels v=33, 45, 57, 93, 127 belonging to the inner well . In the chosen example the laser 
pulse is a Gaussian chirped pulse of energy -Epuise, centered at time tp, with a frequency 

ij{t) = UL+x-{t-tp) = j^i^Lt + m) , (2) 

that varies linearly around the carrier frequency ujl- In Eq. ([2]), we have introduced the 
phase (f){t) of the electric field. The laser is red-detuned from the atomic D2 line at tUat by 
= /l(c<Jat — ^^l), X is the hnear chirp rate in the time domain. The spectral bandwidth 5ui, 
defined as the FWHM of the intensity profile is related to the duration tl of the transform 
limited pulse with the same bandwidth by 5uj = 41n2/r/^ ^ 14.7 cm~^/ri[in ps]. The 
instantaneous intensity of the pulse involves a Gaussian envelope, 
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with a FWHM equal to tq (> tl), the pulse duration [4]. These parameters are related by 
X^T^ = (41n2)^[(rc/ri)^ — 1]. For this pulse, 98% of the energy is delivered in the time 
window [tp — Tc,tp + tc] over the illuminated area a js]. During this time, the instanta- 
neous laser frequency is resonant with the excited levels with a binding energy in the range 
[Al — h\x\Tc, + /i|x|tc], which defines a photoassociation window in the energy domain. 
Assuming the reflexion principle to be valid, this window can be translated into a window 



in the R— domain, considering the outer classical turning points of the set of levels. As 
discussed in [5], for large enough values of the detuning, the photoassociation window in 
the i?— domain is located in a region where the nodal structure of this initial wavefunction 
is energy-independent in the range defined by the thermal distribution, so that thermal av- 
eraging is straightforward. When this is not the case, numerical averaging is possible j^. 
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FIG. 1: (color online)The potential curves for the lower triplet state and the double well 0^(P)3/2 
excited state. Also indicated are the photoassociation windows in the time t (blue), energy A (green) 
and crossing distance Rcii) (red) domains in case of a pulse with negative chirp. 



B. The two-channel coupled equations 

The vibration dynamics in the ground and the excited electronic states is described by 
the time-dependent Schrodinger equation 

H*(t) = (H^oi + W(t))*(t) = ^^^*(t). (4) 
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where ^(t) is a two- component wavefunction describing the relative motion of the nuclei 
in the two channels. The molecular Hamiltonian Hmoi = T + Vei is the sum of the kinetic 
energy operator T and electronic potential energy operator Vei , with components V^round 
and Vcxc on the ground and excited surface respectively . The coupling term is written in 
the dipole approximation: 

W = -D{R) ■ e2E(t) (5) 
£{t) = £ocos{ujLt + ^(t)) 

involving the transition dipole moment of the dimer D{R) and the electric field defined by 
a polarization vector el assumed to be constant and by an amplitude S{t). 

The explicit temporal dependence of the Hamiltonian H is eliminated in the framework 
of the rotating wave approximation. In the present paper, this approximation considers 
the central laser frequency ul, multiplying the radial wavefunction for the nuclear motion 
in the ground and the excited states by exp(— ^^01^^/2) and exp{+iujLt/2) respectively, and 
neglecting the high frequency component in the coupling term. This allows to write the 
radial coupled equations as 

^ I T + V,{R) iy£exp(z0)\ / ^g{R,t) 
\^W^iexp(-z0) T + Ve{R) ) \'^eiR,t) 

where the potentials are now crossing 

VgiR) = V^ound + hLUL/2, Ve{R) = Kxc " (7) 

and the coupling term is defined as 

WL{R,t)exp{-z(f){t)) = -^D{R) ■ e-iSof{t)exp{-'-x{t - tpf) (8) 
In the following calculations, the /^-dependence of the dipole moment was not considered. 




C. Results of the time-dependent calculations 



Time-dependent calculations were performed by integrating Eq.Q, the numerical method 
have been described previously 0, 4]. Due to the low collision energy considered only s— wave 
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scattering has to be accounted. The initial state is represented by a stationary coUision 
wave function, ^g{R,tinit) as illustrated in the upper panel of Fig. [21 corresponding to a 
collision energy E/kB=5A yuK. For a scattering length a ~ 300 ao, the last node in the 
inner region (not considering the pure sine behaviour of the wavefunction for a non zero 
energy), hereafter referred-to as "last node" is located at Ra = 310 ao(~ a), and the last 
but one at -Riasti= 82.5 ao. In an energy range close to threshold, all nodes with R < -Riasti 
are energy-independent. The cold atoms are illuminated by a chirped laser pulse , referred 
to as in j^gfs 

, |6||, with central detuning = 0.675 cm ^, duration (FWHM) 
Tc = llOps, stretching factor fp = tc/tl= 1.91 associated with a transform limited pulse of 
duration r/^=57.5ps and maximum intensity Il= 120 kW cm~^. The central frequency Ul is 
resonant with the level f =153, for which the outer turning point is located at Rl = 148.5 ag 
( corresponding to a maximum of the wavefunction, beween -Riasti and -Ra)- For a negative 
chirp, the photoassociation window spans the levels v= fmax = 159 till Vjam =149, bound by 
0.456 to 0.869 cm^^, with a classical vibrational period Tvib from 1095 till 635 ps. In the 
i?— domain, the window extends from i?min = 135 till -Rmax = 176 ao. 

The conditions for an impulsive approach are that the pulse duration is shorter than the 
vibrational periods (see Ref. j^). At the intensity considered, during the pulse many levels 
are populated, but after the pulse only the levels between Vmm and fmax remain populated 

0. 



D. Momentum kick on the ground electronic state 

The generation of the depletion hole on the ground a^S^ electronic state can be under- 
stood by employing a phase amplitude representation both for the ground and excited state 
functions: 

^,AR,t) = A,,,(i?,t)exp[^^^M] = A,,,(i?,t)exp[^£^^^^^l^], (9) 

where Pg{R,t) can be interpreted as the local momentum in the ground state [3, Q]. The 
evolution of the amphtude Ag{R,t), displayed in Fig. [21 suggests the following remarks: 

• Before the pulse, we see in the upper panel the initial state stationary collision wave- 
function |\l/c,(i?, tinit)| described above, and computed as a unity normalized eigenstate 
in a box of length L=19 250 slq. In order to optimize the population transfer, we have 
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FIG. 2: (color online) Upper panel (blue line): amplitude Ag{R,t = Unit) of the initial scattering 
stationary wavefunction for a pair of ground state cesium atoms colliding with energy E = ksT, T = 
54iJ,K, and zero angular momentum, in the lack of any electromagnetic field. Note the position of the 
"last node" at Ra and the last but one at Riasti- Upper panel(red line) : amplitude Ag{R,t = tp) 
at the maximum of the pulse t = tp : the pulse has carved out a hole in the initial wavefunction. 
This hole is indeed located in the region [i?mirn -Rmax] of the photoassociation window, delimited 
by the two vertical lines. Next panels : after the pulse, the ^g{R,t) wavepacket moves inwards : 
note the increase of the maximum value of the amplitude in the inner region. At t — tp ~ 850 ps, 
the wavepacket has been partly reflected by the ^T,^ inner wall located at R ^ lOoo. Note that 
in the excited state the classical vibrational half-period for the photoassociated levels is in the ~ 
[300,600] ps range 
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chosen a detuning such that the photoassociation window [135 , 176 ao] is located in 
the region of the "last" maximum of \'$g{R,t = 0)|. 

The laser pulse, maximum at t = tp, carves out a depletion hole in the initial 
wavepacket. The location of this hole is consistent with transfer of population to 
the levels f =149 to 159, with outer turning points from 135 to 176 ao, located within 
the photoassociation window. Note that since J^^^ \'^g(R,t = 0)\'^dR = 0.0025, only a 
very small fraction of the probability density in the initial state is transferred to the 
excited state. This very small value is due to the fact that we use a very large box. 
Calculations using energy- normalized initial state have been reported elsewhere js, 4|, 
as well as determination of the absolute value of photoassociated molecules 

After the pulse, this hole starts moving to smaller internuclear distances, due to a 
negative "momentum kick" (see Ref. and Sec JIIII below). Whereas the motion is 

taking place in a region where both Vg and ^ are negligible, the classical velocity is 
quite impressive : for instance, the new maximum created on the left side of the hole 
and initially located at -R=110 ao moves to -R=90 ao with a velocity ~ 4.2 m s^^, two 
orders of magnitude larger that the classical velocity of the colliding atoms. 

The motion is restricted to a localized part of the wavepacket : during the time-period 
considered here, the position of the "last" node at Ra= 314 ao, and the wavefunction 
at distances larger than Ra, are not modified. 

In contrast, we observe an increase of the probability density at distances shorter 
than the hole position, hereafter referred to as "compression" of the wavepacket. The 
amplitude of the secondary maximum created on the left side of the hole increases and 
has grown by a factor of 2 when it reaches R ~90 ao. Therefore, at times t—tp >350 ps, 
(the optimal value being 850 ps ), the wavepacket is well adapted to photoassociation 
with a second pulse, red-detuned relative to the first one, populating in the excited 
state vibrational levels with an outer turning point around 90 ao. 

At t — tp=850 ps, the wave packet is reflected by the inner wall of the a^E+potential, 
and starts moving outwards. Note that in the excited state the classical vibrational 
half-period varies from 318 for Vmin to 548 ps for Vmax- 
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We may therefore conclude that after the pulse the motion of the wavepacket in the 
ground state is governed by a timescale similar to the vibrational motion in the pho- 
toassociated levels of the excited state. 



We next display in Fig. [3] several snapshots for the variation of the phase derivative 
^^J^^ = Pg{R-,t)/h (see Eq. OH])), where pg(R,t) can be interpreted as the semi classical 
momentum [tI. Just after the pulse, significant momentum is created in the region of the 
photoassociation window, where the potential and its gradient are negligible. The computed 
value ~ 0.20 (ao)'^ corresponds to a classical velocity of 3.6 m s~^, same order of magnitude 
as the very rough estimation of 4.2 m s~^ given previously. As expected, a strong momentum 
increase is observed in the short range region ( i? < 30 ao), where the potential is strongly 
attractive and the wavefunction rapidly oscillating. This will be analyzed further in Sec. 

m 



E. Advantage of the compression effect for photoassociation with a second pulse 

To explore further the possibilities offered by the compression of the wavefunction we have 
represented in Fig. Hlthe time- variation of the Franck-Condon overlap between t) and 

the wavefunctions of the bound vibrational levels in the excited state. Almost all levels in 
the external well of the 0~{6S + 6P3/2) excited potential curve can be efficiently populated. 
The deepest level populated is v=29, v^xt = 4 with a binding energy of 70 cm~^, an inner 
turning point at i?=30 ao, and a good Franck Condon overlap with the v"= 42 and 33 levels 
of the a^S^ potential, respectively bound by 3.7 and 20.5 cm~^. Population of such levels 
might be interesting for the implementation of the stabilization step. 



F. Redistribution of population in the ground state : formation of pairs of hot 
atoms and of halo molecules 

As a result of the short pulse, the population in the ground state is redistributed. Indeed, 
the photoassociation reaction ([T]) is accompanied by a redistribution process, 

2Cs(6^, F = A){E = keT) + h{uj{t), ) ^ 2Cs(65, F = A){E' = keT') (10) 



9 



n dSg(R,t)/dR 

0.25 
-0.25 



-J PA window 




FIG. 3: (Color online) Same as Fig. [2] for the variation of the semiclassical momentum pg{R,t) = 
dSg{R,t)/dR, where Sg{R,t)/h is the phase of the wavefunction in the ground state (see Eq. 
([9]) in text) . Just after the pulse, the value ~ 0.20 (ao)^^ of the momentum created in the 
photoassociation window is consistent with a classical velocity of 3.6 m s~^. As expected the 
momentum is increased in the short range region where the potential Vg{R) is not negligible. 

producing pairs of ground state atoms with different energies ( in particular, pairs of "hot" 
atoms that may leave the trap) and by an association reaction. 



2Cs(65, F = 4) + h{uj{t)) Cs2(^S+(65 + 63); v", J) 



(11) 



with formation of bound molecules in the last vibrational levels of the a^S+potential. 

The redistribution into the stationary levels $g(f") of the ground state potential is an- 
alyzed in Fig. [5] as a function of the final collision energy {E' = ksT') by projecting 
the wavepacket \l/g(i?, t) after the pulse on eigenstates (energy-normalized continuum lev- 
els and bound levels) of the (T + Vg{R)) Hamiltonian. The calculations are considering 
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FIG. 4: (Color on line)Interest of the compression effect. Intermediate panel : wavefunction 
t = tinit)| of the initial collision state and ground state wavepacket \'^g{R, topt = ip + 850) |, 
850 ps after the pulse maximum, when the compression effect is maximum. Upper panel: integrated 
population \^g{R, showing the increase of the probability density at short distances. Lower 
panel : variation of the overlap | < ^e{v)\'^g{R,t) > p with the eigenfunctions for all the bound 
vibrational levels v in the 0~(65+6P3/2) excited potential curve, as a function of the binding energy 
Eb{v), for t = tinit (black) and t = topt- Several dips visible in the overlap should be attributed to 
levels in the inner well ( v=33, 45, 57, 73,93, 127), which are indicated in red. 

a unity-normalized initial state. We have studied the sensitivity to the pulse parameters. 
The redistribution is independent of the sign of the chirp; in contrast, when the maximum 
intensity of the pulse II (see Eq. (|3])) is increased, the width of the redistribution function 
is clearly reduced. Note that the area below the curve is roughly conserved, changing from 
0.0097 to 0.011 when the intensity increases by a factor of 9. 

Population is also transferred efficiently to bound levels of the ground triplet state, pri- 
marily to the last bound level v"=53, creating a halo molecule. This effect is similar to the 
population of the last bound level by sweeping an optical Feshbach resonance. The num- 
ber of ground state halo molecules (0.004) is the same order of magnitude as the number 
of photoassociated molecules (0.002) as can be seen by comparing the left and right lower 



panels of Fig. [5], in accordance with Ref. js]. It should be noticed that the population of 

11 



d N3,*j„ /dE = i< i-j <i>^(E') >r 




100 2 

N„-^E_,= l<1'„J*,Ji/) 
PA window 



em 



0.005 -NVj^,,= |<>I'J*g(0>| 





°i — i — 

FIG. 5: (Color on line)Population redistribution in the ground state : analysis of the wavepacket 
after illuminating with the pulse P}_^^, with negative chirp (black circles), and P^'^ (red squares). 
Redistribution after a similar pulse with higher intensity, changing the maximum from II to 91 l{ 
green up triangles, x < 0) blue down triangles, x > 0)- Upper panel : redistribution in the 
continuum, as a function of the energy, for an initial state -E'j„jt=54 ^ K. The continuum levels 
are energy normalized. The results do not depend upon the sign of the chirp. Right lower panel 
: population transfer to the last bound levels of the ground triplet state, as a function of the 
vibrational index : only the last level is populated. Left lower panel: for comparison, population 
transfer to the bound levels of the excited state, as a function of their binding energy. 



the photoassociated molecules is spread over more than 30 vibrational levels. 

When the sign of the chirp is changed to a positive value, the population of the last bound 
level and the redistribution into continuum levels is not affected. In contrast, the photoas- 
sociation efficiency is different (see left panel of Fig. For x < 0? the levels populated in 
0^(6S+6P3/2)are within the photoassociation window, and the population depends weakly 
upon II. For x > 0, the distribution extends to levels outside the photoassociation window, 



and this effect increases with intensity. As discussed in Ref. 1^, 11], for x < the instanta- 
neous frequency "follows" the motion of the wavepacket, so that the population is brought 
back to the ground state and up again. This phenomenon is labelled "multiple interaction" 
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in Ref. For x > O5 this recycling effect is much wealcer : once a level is populated, the 
instantaneous frequency is no longer resonant with its energy, and the population remains 
in this photoassociated level. This recycling effect will be discussed further in the following 
sections. 



III. THE CURRENT DENSITY VECTOR AS A TOOL FOR INTERPRETATION 



A. Local current density vector and local momentum 

In order to get better insight into the hole dynamics one has to follow both the population 



For the 



depletion as well as the local value of the momentum induced by the laser field 
problem of photoassociation, a global quantum analysis of the total changes in population 
and in momentum due to the pulse has been discussed recently in Ref. 91]. This formalism, 
based on Heisenberg equations, is adapted to the implementation of local control theory. 
In the present work we use a local framework closer in spirit to the phase amplitude sepa- 
ration Eq. (ini) and based on the change in time of the radial component of the probabilty 
current density P, 18|. The analysis is carried out for both the ground and excited state 
components of the wavefunction as defined in Eq. ([6]), 

3g/e{K, t) - ^i^9/el^. t) "^g/eKii, t)\. [12) 

The local momentuni, defined as semiclassical momentum in Eq. ([9]), and called "state 
momentum" in Ref. [lO(], is related to the probability current density 

where the phase and amplitude have been defined in Eq. ([9]). 



The formulation is related to the classical limit of Bohmian dynamics |12!], where one 
considers on each channel the motion of a fiuid of non interacting particles with probability 
density |^3,eP, probability current density jg^e{R,t), driven by a force Fg^e{R) — '^^^'^ — 



dt 



--dt 



+ -^^)pg^e{Rit)- This aspect will be discussed in a forthcoming paper [l 
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The time derivative of the two quantities mjg{R,t) and mje{R,t) reads 

^djgi^ = Fg{R)\Ag{R, t) ^ + lCg{R^ t) + C{R, t) (14) 

^djA^ = F,{R)\A,{R, t)\' + IC,{R, t) - C{R, t), (15) 
ot 

where m is the reduced mass, and we have defined 

• the classical force for each potential 

F^/eiR) = (16) 

so that the first term in Eqs. ( iMllTSi) depends both upon the classical force and the 
local probability density y4^yg(i?, t) in the channel considered, defining a density of 
force Fg/e{R)Al/^{R,t) . 

• a kinetic term due to the non-hermiticity of the local kinetic energy operator 

This complicated term, arising from the fact that the time- variation of j{R,t) results 
both from the time- variation of the local momentum p{R, t) and of the local probability 
density \A{R,t)\'^, will be discussed further in a forthcoming paper [isl . 

• a third term due to the coupling by the laser field, which cancels when the laser is off, 
and has an opposite sign in the ground and in the excited state. 

C{R, t) = ^mR, t)^{WL exp(-z0(t))v^,(i?, m-^[Wl exp (^0(t)) ^;{R, t)^^^] 

(18) 

It depends upon the coupling (amplitude Wl, phase 0, and therefore sign of the chirp) 
and also upon the phase difference between the wavepackets in the ground and in the 
excited state. Introducing the momentum operator P = —ihd/dR , one may rewrite 

CiR, t) = -1 [S(vl/*(i?, t)PiWL{t) exp(-#(t))^,(i?, t))-3(vl/;(i?, t){Wlit) exp(+z0(t))Pvl/,(i?, t)), 

(19) 

which is close to the notation of Ref. where it was termed the momentum exchange. 
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Note that the density of force C{R,t) is zero if the laser field is zero, of course, but also 
when the local population either in the ground or the excited state is 0. 



Another set of equations, complementary to Eqs. fll4lll5p describes the variation of the 
local probability density 

dAg{R,ty _ djgiR,t) 2^ 



m- 



dt OR h 

dA,{R,tf dUR,t) , 2 



-^{^:{R,t)WL{t)exp{-t<p{t))^g{R,t)), (20) 



-%^:{R, t)WL{t) exp{-t<p{t))^,{R, t)). (21) 



dt OR h 

When the laser is off, Eqs. fl20|21l) are the relations of flux conservation on each channel. 
During the pulse, transfer of population between the two channels is taking place. 



B. Integrated and partially integrated values 

With the concept of the photoassociation window in mind we can differentiate between 
the local and global influence of the pulse. For the global effect, we have considered the inte- 
grated values, over the total box of length L. We define for the probabiltiy current densities 
jg/e{R,t), the position probability densities \Ag/e{R,t)\'^, and for their time derivatives, the 
integrated functions: 

mlg/e{t) = ! mjg/e{R,t)dR (22) 
Jo 

m — = m I — dR (23 

dt Jo dt ^ ^ 

Ng/e{t) = j \Ag/,{R,t)\''dR (24) 

Jo dt ^ ^ 

mlg/e{t) is the current and Ng/e{t) is the number of particles. In Ref.j^, similar quantities 
are obtained directly from the Heisenberg equations. The quantities corresponding to the 
integrated value of )Cg^e{R,t) then disappear. In fact, in the present approach, when local 
time-dependent quantities are integrated over the whole box, the hermiticity of the kinetic 
energy operator results into JCg e{R,t)dR = and one obtains 



m 




dt 



/ Fg/,iR)Al/^iR,t)dR± / C{R,t)dR (26) 
Jo Jo 
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so that 

m I ^(Wi) + Ie{R,t) ^^ ^ {F,{R)\A,{R,t)\' + F,iR)\MR,t)\')dR (27) 

where the term due to laser couphng has disappeared. 

However, for the analysis of the photoassociation effect, where the process is taking place 
within a photoassociation window, a local view is consistent with partially integrated values. 
We have therefore considered also the time dependent quantities Igjl^it) and N^/l^{t) 
defined as integrals 

fR2 



J,/e{R,t)dR (28) 
N^Jeit)\t = / Al/eiR^t)dR (29) 



taken on a limited range of internuclear distances [Ri,R2], that could be [Rmin, Rmax]- In 

rR2 



the general case, the influence of the kinetic term [j^^ ICg^e{R,t) will manifest itself. This 



will be discussed in detail in a forthcoming paper [1 

In the present work the limited range was chosen as [0,-Ra], so that j^"^ K,g^e{R,t) is 
negligible. Ra is the position of the "last" node, which stays fixed during and after the V^'^ 
photoassociation pulse (see Fig. [2]), as discussed above in Section [Tl CI The justification for 
such a choice is that at large distance, where the potentials are negligible and the dipole 
moment Z)(i?) constant, the wavefunction in the excited state merely reflects the initial 
wavefunction 3|, l^] ■ Writing 

^g/e{R,t) ^ agle{t)^!g{R,t 

init j 

where ag/e{t) is a complex number, such that ae(t ^ tc) = 0, it is clear that there is no 
i?-dependence of the phase Sg/e{R)- The analysis of the photoassociation dynamics can 
then safely be restricted to the R < Ra range. 

In the following we shall simplify the notations by writing 

C(^) = I^Jeit) 1^ (30) 
K/eit)=^7e\t)\0^ (31) 
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IV. ANALYSIS OF THE RESULTS 



A. Analysis of the results for the pulse V^'^ 

A comparison of the time- variation of mlg/e{t), l^g/e{t), and their time derivatives, as 
well as of the partially integrated quantities ml^^'^^{t), N^Jl^{t) are reported in Fig|6l We 
have also defined a "partially integrated" force, 

F7f (0, Ra) = r Fg/,{R)Al/,{R, t)dR. (32) 

J 

A look at Fig. [6] leads to the following remarks: 

• During the time window tp — tc < t < tp — tc, there is a large population exchange 
between the two channels : in the lower figure of the left pannel, '^^^^^^ = — ^^^^ 
displays a variation roughly identical to the shape of the pulse, apart from small 
oscillations. The range of distances i? > 314 contributes very much to this population 
transfer, and we see that the variation of the partial population, ^^^^^ is much 
smaller than '^^^^^^ ■ 

• After the pulse, the population transfer is restricted to the photoassociation window, 
as displayed below in Fig. IIV A[ 



In contrast, there is no momentum exchange at large distances, and we find Ig/e{t) 

^par 

9/e 



I^Teit) ■ Indeed, in the region R > Ra the phase is i?— independent on both 



channels(^^2|m = q). 

At the very beginning of the pulse, for times t < tp—Tc, when, due to low laser intensity 
the (small) population transfer is not adiabatic, the ground state wavepacket gains a 
positive momentum . The total momentum Ig + Jg is conserved, the contribution of 
the force being negligible [l|. 

During the time window [tp — Tc,tp + tc], there is an exchange of probability current 
between the ground and the excited state. The two currents Ig and Jg are oscillating 
with opposite phase, in agreement with Eqs.( l26l) . Therefore the wavepacket created 
in the excited state moves locally back and forth during the pulse. 
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FIG. 6: (Color on line) Photoassociation with the pulse "Pi^^, with negative chirp and small 
detuning, maximum at t = tp=350 ps: the time-window [tp — Tc,tp + tq] is indicated by the 
horizontal arrow. Left column : in the upper figure is represented the integrated value of m 
times the probability current density mIg{R^t) > (red line) in the ground and in the excited 
state mIe{R,t) > (black line), and of their sum as a function of time (green line). Note the Rabi 
oscillations during the pulse. A small current in the ground state at the beginning of the pulse 
is drawn as enhanced by factor 100. In the lower figure, the time variation of the population is 
displayed: the calculations have been done by summing over the full range of distances (black lines) 
and by restricting to R < Ra =~ 314 oq (red lines). It is clear that a large amont of population is 
transferred at large distances during the pulse. In the right column, the time variation nT-—^^^^ (see 
Eq. [26] in text) in the excited (upper figure) and in the ground state is compared to the 'classical' 

^par 
9/e 



term F''^'^*(0, i?a)i showing a strong enhancement during the pulse. Note that the classical force is 



negligible in the ground state but not in the excited state. After the pulse both curves follow the 
classical force 
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The sum of the probabihty current on the two channels is not oscillating, in agreement 
with Eq. fl271) . Due to the large acceleration in the excited potential (see the upper 
figure in the right pannel), this sum current is smoothly decreasing. 

For tp — Tc < t < tp the two channels are equally sharing the classical acceleration 
due to the asymptotic behaviour of the 0~(6S+6P3/2)potential, so that Ig{t) and 
Ie(t) are oscillating with opposite phase below and above their mean value. 

For tp < t < tp + Tc the current density gain in the ground state becomes larger than 
in the excited state : the acceleration is no longer equally shared between the two 
channels. 

During the pulse, the time derivatives m^^^^^^ are much larger than the classical forces 
Fgie{R)A^^I^{R, t)dR in Eqs. ( j26l) . The photoassociation dynamics is dominated by 
the laser-induced force arising from the term J^^" t)dR. 

After the pulse, the inner part of the wavepacket in the excited state is acceler- 
ated towards short distances, due the classical acceleration. The pulse has 
been optimized for focussing this vibrational wave packet half a vibrational period 
after the pulse maximum [3| : for the levels considered, 318 ps < Tmb{v)/2 < 550 
ps. The refiexion on the inner wall occurs during a very short time, in the range 
tp-rc + T^ib{vmax)/'^,tp + Tc + T^ib{vmin) f^, ^-6. [778,790] ps. In Fig. Owe may see 
a sudden change of the sign of Jg at {t = t\ 

oc ~ 784 ps. 

In contrast, as can be seen in the lower figure of the right pannel, the classical accel- 
eration is zero in the ground state. After the pulse, Ig remains constant during a long 
time delay of ~600 ps. Since the population remains constant, this corresponds to a 
constant velocity. The part of the wavepacket that is located dX R < Ra is moving 
to short distances with constant velocity for a duration of ~600 ps, in a region where 
there is no potential gradient. We see then the evidence for a negative momentum 
kick due to the laser pulse. 

For t — tp > rc+600 ps, the wavepacket in the ground state reaches distances where 
the a^S+potential is no longer negligible : it is accelerated at ti 1000 ps, then partly 
reflected by the inner wall at t ~ 1100 ps. At t = topt = 1300 ps {t — tp= 950 ps), the 
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mean velocity vanishes, which corresponds to the optimal time for reexcitation by a 
second pulse ( see Fig. H]). 

The oscillations which appear during the pulse in the time variation of the populations or 
of the current are analyzed in more detail in Fig. IIV A[ These Rabi oscillations are induced 
by the laser field : indeed, for the V}^"^ pulse the Rabi frequency at resonance, evaluated as 
n{t) = \\WL\f{t)) satisfies 

r tp+TC 

/ n{Rc,t)dt 5.3U, 

J tp-TC 

, which is the same order of magnitude as the number of oscillations visible in the figure. 



• The n phase difference between the oscillations of the populations Ng and results 
from the conservation of population . 

• In a given channel, the 11 phase difference between the oscillations of the proba- 
bihty current Ig and /g results from the ± jC{R, t)dR laser couphng term in Eq. ( !26l) . 



• In each channel, the 11 phase difference between the time variation of / and points 
out that both in Eqs. fll4lll5p and in Eqs. fl20|2ip ,the laser coupling term dominates, so 
that the change in the current is mainly caused by the change of population. In other 
words, the exchange of population between the two channels induces the momentum 
kick in the ground state: the population transferred back to the ground state causes 
a negative momentum kick. It is the Rabi cycling phenomenon, and the efficiency in 
recycling population that seems the key factor in inducing a motion of the wavepacket 
towards short distances. 

The average value of the momentum gained during the pulse by the inner part of the 
ground state wavepacket can be evaluated from 

Nr'(t) 

This quantity is represented in the left upper pannel of [HI The Rabi oscillations observed 
in the time- variation of I^^^^{t) and NP'^''*(t) being of opposite phase, they disappear almost 
completely in < PP'*''*(t) >, which proves again that the oscillations in I^^^*'{t) are almost 
entirely due to the population. The mean negative velocity gained by the ground state 
wavepacket is ~ 2.7 m s~^, in agreement with the previously estimated value. 

20 



< >= J^J,', (33) 
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FIG. 7: Evidence for a momentum kick in the ground state. The laser pulse is P^^^ maximum at 
time t =350 ps. Top figures : Time- variation of the probability current Ie(t) = Ie^'^{t) (left) and 
Ig{t) = Ig^'^^{t) = (right) on the two channels . Middle pannel : time variation of the partial popu- 
lations N^'*'^*(t) = and Ne'*'^*(t) = in the range [0,R = 314] of internuclear distances. Lower pannels 
: time-derivatives of the same quantities. The vertical lines are drawn to emphasize the phase of 
the Rabi oscillations for each quantity. The pulse is present during the time window [240,460]ps. 
After the pulse, whereas the acceleration in the ground state is negligible {dlg[t) /dt=0), the average 
momentum is negative and constant, manifesting the momentum kick in the ground state. 



B. Sensitivity to the sign of the chirp parameter and to the intensity 

The sensitivity to the sign of the chirp is illustrated in the right pannel of Fig. [8], where 
are reported calculations done for the photoassociation pulse similar to the previous 
one but for x > 0. The maximum population transferred to the excited state during the 
pulse is much smaller than for a negative chirp, the amplitude of the oscillations being also 
smaller. However, after the pulse, the population remaining in the excited state is 0.003, 
larger by a factor 1.5 than for x < 0. Looking at the mean value of the population, we 
see that recycling occurs for the pulse whereas for a positive chirp a regular increase 



of the popu 
interactions 



ation is observed. This recycling effect can be attributed to the multiple 



10 



111 ] occurring for a negative chirp, when the motion of the two wavepackets 
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time window 





1000 
t(ps) 

FIG. 8: Average momentum and influence of the sign of the chirp . Left pannel: Photoassociation 
with the pulse 7^1^^, with negative chirp : Comparison between expectation values of the momentum 
< Pg^^^{t) > and the probability current Ig{t) ,the population in the ground state Ng'^'^*(t) in the 
region R j 314 ao(red lines) . The quantities for the excited state, < Pi^^^{t) > , Ie{t) and Ng'^'^*(i) 
are drawn with black lines. . Right pannel: same quantities for photoassociation with the pulse 
PP^, with positive chirp. 



is "following" the pulse, and not for a positive chirp. 

After the V^"^ pulse, the wavepacket in the excited state is spreading, since there is no 
focussing effect : the reflexion on the inner wall takes place during a longer time interval 
than for x < 0, from tp - tq + \T^ib{vmin)= 560 ps to tp - tc + \T^ib{vmax) = ^^^^ ps. 
Starting from t= 800 ps, the large spead of the wavepacket gives rise to interferences 
between components moving in the opposite direction. 

The mean momentum gained by the inner part of the wavepacket after the V^"^ pulse 
corresponds to a mean velocity of 2.1 m s^^, a factor 1.3 smaller than in the case of the 
negative chirp. 

It seems that the momentum kick is favoured by important Rabi oscillations and a smaller 
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photoassociation efficiency. 



The effect of the laser intensity is discussed in Fig. O When the intensity is increased by 
a factor of 9, the number of Rabi oscillations during the pulse increases by a factor of 3 as 
expected. The transfer of population after the pulse NP'^''*(t > t + tc) increases from 0.0024 
for X < to 0.0038 for x > 0. The population in the excited state is now only a factor of 
1.15 larger for x > 0. : this can be explained by the occurence, at larger intensities, of a 
recycling effect not only with negative but also with positive chirp. 

When considering the larger intensity and a negative chirp, the momentum kick in the inner 
part of the ground state wavepacket is increased by a factor of ~ 2. The inner part of the 
wavepacket moves to short distances with a mean velocity of ~6 ms~^ ( compared to 2.7 
previously) : the reflexion on the inner wall of the a^E^ potential, with /p=0, occurs ~ 400 
ps after the pulse maximum, to be compared with the value 950 ps in the previous case. 
For positive chirp, the momentum kick in the ground state is now very close to the values 
obtained for x < 0. However, since in the latter case the population NP'^''*(t > t+Tc) is larger 
by a factor of 1.5 than for x > 0, a more important compression effect is to be expected for 
negative chirp. 

The previous analysis shows that the relative gain in momentum on the two channels can 
be modified by changing the parameters of the pulse. The way of optimizing the latter is 
an open problem. 



C. Local and non-local effects 



During the pulse, the variation of the local probability current density on each channel, as 
described in Eqs. fll4|15p . as well as the probability amplitudes (see Eqs. ( I20|2T1) depend upon 
coupling terms involving the laser coupling WL{t) exp(— , and either 
or t) '^g{R,t). The important quantity is then the phase difference, 

Sg{R,t) S,{R,t) 



h h 



m (34) 



which depends both upon t and R. In our problem, the initial wavefunction t) is 

delocalized, and extends over the full range of internuclear distances. Therefore during the 
pulse the excited wavepacket \l/e(-R, t) will also be delocalized. However, for very large dis- 
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FIG. 9: (Color on line)Effect of the sign of the chirp and of the intensity of the photoassociation 
pulse. Time-variationof the current mig in the ground state (upper figures) and m/g in the excited 
state (lower figures) . Left column, calculations with the pulses Pl^^ (black lines), 'P^'^^(red lines) 
as in FiglHl Right column : same calculations, increasing the coupling in Eqs.® by a factor of 3, 
or the intensity II in Eq.([3|) by a factor of 9 and considering negative (black lines) and positive 
(red lines) chirp. The increase of the momentum kick in the ground state is clearly manifested 



tances, neither Sg{R,t) nor Se{R,t) present a R dependence. In Fig. [TOl the i?-variation of 
^^^^hdR^'^'' represented for various times in the vicinity of the pulse maximum. For distances 
up to 250 ao, very rapid changes of the phase occur, and a very complicated /^-dependence 
is manifested, manifested. This suggests that the choice of /^-integrated physical quantities, 
as used in Refs. jo], 1^, to optimize the pulse, may not be a very efficient procedure when 
the delocalized character of the initial wavefunction is introduced in the model. 
In the "single pulse Raman-like walk scheme" of Ref . [lo| . the initial state is a localized gaus- 
sian wavepacket, and the control parameter is the laser instantaneous frequency, designed 
to match the collision dynamics at distances R <40 ao, by inducing multiple coherent up- 



24 



h"'dS„(R,t)/dR 




tp+100 



V150 



lastl ■ 'a ■ "lastl 

FIG. 10: (Color on line) Snapshots for various times during the pulse of variation of the local 
momentum see Eq.Q as a function of the internuclear distance R. This figure is 

giving more detailed information than FiglSl the strong variation of the phase as a function both 
of the distance R and of the time t is demonstrated. The position Ra of the "last node" in the 
initial wavefunction, and Riasti of the previous one are indicated by arrows; 



ward and downward transitions. It seems that for this "matched spectrum pulse", a lot of 
population is redistributed into continuum levels of the ground state, creating pairs of hot 
atoms rather than the target bound levels in the ground state. In Ref. js], the global phase 
of matrix elements such as < '^e{R,'t)\D{R)\'^g{R,t) > and < ^e(-R, ^)|-D(-R)| 



dR 



> is 



the control parameter. The phase of the field is then following the variation in the phase of 
the expectation value to be controlled. It then seems that for pulse duration larger than 1 
ps , the efficiency of the control is dropping. 

The choice of a spatially global operator such as the total population or the total momentum 
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can overlook the finer details of controling a localized quantity. This is equivalent to the 
i?-integrated quantities in the present paper, which loose a large part of the phase informa- 
tion. Therefore, the partially integrated quantities (I29f3ip are a better choice, enabling a 
better optimization of the compression effect. In the next paper the current approach 
will be linked to a different set of spatially semilocal operators. 

However, from the present analysis of the partially integrated quantities, we have drawn a 
few conclusions that might help to optimize the compression effect. 

V. CONTROL OF THE COMPRESSION EFFECT 

From the discussion of Section IIVI it seems that the creation of an important flux of 
population to short distances is favoured both by an important population cycling during 
the pulse and by a good transfer of this population back to the ground state. In other words, 
the optimal pulse should be as close as possible to a {2n + 1)11 pulse, with a negative chirp, 
and a large number n of cycles. This can be achieved with a pulse that is "following" the 
excited wave packet. We have tentatively changed the parameters of the photoassociation 
pulse to check whether it is possible to optimize the compression effect. Starting from a 
transform limited pulse with the same spectral width as previously ( 5ct;=.255 cm^^, Ti= 
57.5 ps) and with the same central frequency, the best results were obtained by increasing 
the linear chirp parameter x ^^id hence the temporal width. The optimal value was tq = 
376.13 ps, instead of 110 ps in the previous calculations. The duration of the pulse then 
becomes comparable with half the vibration period (T^jfc/2)in the excited state. In order 
to make results comparable, the eo factor in [8] was mutiplied by a factor compensating 
the variation of a/ tl/tc factor in f{t) in order to keep constant the maximum value of 
the coupling. The energy Epuise is then multiplied by 6.25. The new pulse is closer to a 
(2n + l)n one, the number of Rabi oscillations reaches up to 24, the population transferred 
to the excited state has markedly decreased (3 10~^ instead of 2 10~^). 
The compression effect is now much more important, occurs earlier and is optimum at 
iopt = tp+ 350 ps. The wavepacket at this time is displayed in Fig{TT|, demonstrating a 
spectacular increase of the compression effect as compared to the previous results, and 
suggesting that the direction is the good one. 
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FIG. 11: Improvement of the compression effect presented in Fig. HI using the new pulse described 
in Section |Vl with a parameter tq comparable to the half vibration period in the excited state. 
Intermediate panel : wavefunction \^g{R,t = tinit)| of the initial collision state and ground state 
wavepacket \'^g{R,topt\ when the compression effect is maximum, after illuminating with the pulse 
"Pq^^ (black line,fopt = tp+ 950 ps) and the new pulse (red line, topt = tp+ 350 ps). Upper 
panel: integrated population \"^g{R, t)\'^, showing the increase of the probability density at short 
distances. Lower panel : variation of the overlap | < <I>e(u)|^'g(i?, t) > p with the eigenfunctions 
for all the bound vibrational levels v in the 0~(65 + 6P3/2) excited potential curve, as a function 
of their binding energy Eb{v), for t = tinit(blue) and t = topt (black and red line). The compressed 
wavepacket is even better adapted for further photoassociation. 

The analysis of the results shows that the model of adiabatic transfer within a photoasso- 
ciation window is still valid, but the impulsive approximation is no longer valid; new physical 



effects are present, which will be analyzed in a forthcoming paper 13|. 



VI. CONCLUSION 



The present paper was devoted to the analysis of the formation of a dynamical hole in 
the wavepacket describing the relative motion of a pair of cold atoms in presence of a pho- 
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toassociation pulse, and of the evolution of this wavepacket after the pulse. As a case study, 
photoassociation of Cs into loosely bound levels of the outer well in Cs20~(6S+6P3/2)was 
considered, using a pulse in the picosecond range, with a negative linear chirp parameter, 
designed in previous work {s, 6] to create an excited wavepacket with focussing properties. 
In the initial state, the wavepacket is a stationary continuum wavefunction at the collision 
energy corresponding to 54 fiK. Photoassociation with such a pulse can be interpreted as 
an adiabatic transfer of population within a photoassociation window, while the impulsive 
approximation is valid, since the motion of the nuclei can be neglected during the pulse. 
The numerical calculations show the creation of a dynamical hole in the region of the pho- 
toassociation window. The latter is created at distances where the ground state potential 
is negligible, while the excited state potential, with asymptotic behaviour, is giving to 
the excited wavepacket an acceleration towards short distances. 

After the pulse, the inner part of the ground state wavepacket is shown to move towards 
short distances, at a velocity typical of the vibrational motion in the excited state, due to 
the negative momentum kick gained from interaction with the laser field. This leads to 
a compression of population at short range, so that a second photoassociation pulse, red- 
detuned relative to the first one, and conveniently delayed, could be designed to bring much 
population to deeply bound levels of the excited state. Another signature of the dynamical 
hole is the redistribution of population in the ground state, with formation of pairs of hot 
atoms and population of the last bound level. 

We have analyzed the momentum kick phenomenon by considering the time-evolution of 
the position probability density and of the probability current density. Spatially integrated 
values, giving the population and momentum on each channel, were considered, as well as 
partially integrated values, considering population below a certain distance. This procedure 
is necessary due to the delocalized character of the initial wavepacket. The various terms 
included in the treatment have been discussed in the framework of a phase amplitude for- 
malism: besides the classical force and the force due to the laser, a third term has been 
pointed out. 

The analysis for a pulse in the 100 ps range, widely used in previous work, with negative 
chirp, shows that, during the pulse duration, population and momentum exchange between 
the two channels occur due to Rabi cycling, so that the acceleration from the excited channel 
is transferred to the ground channel. After the pulse, the momentum gained by the inner 
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part of the ground state wavepacket depends upon the number of Rabi cycles, upon the 
sign of the chirp ( it is smaller in case of a positive chirp) and upon the amount of popula- 
tion transferred back to the ground state. It seems to be maximum for a pulse as close as 
possible to a (2n+l) IT pulse, with which no photoassociatcd molecules are formed. In con- 
trast, when considering a similar pulse with positive chirp, the photoassociation efficiency 
increases while the momentum kick is reduced. 

Following these ideas, we have tentatively proposed photoassociation by a longer negatively 
chirped pulse, with duration similar to half the classical vibrational period in the excited 
state, so that the impulsive approximation is no longer valid. Such a pulse is designed to " 
follow" the motion of the excited wavepacket during half the vibrational period. A strong 
increase of the momentum kick and of the compression effect in the ground state was indeed 
obtained, demonstrating that the direction is promising. Further analysis will be given in a 
forthcoming paper. 
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